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Tvfo  aircraft  refueling  syste:;n  configurations  were  sinulated 
using  HYTRAN,  an  existing  transient  flow  analysis  computer  program. 
Transient  pressure  response  subsequent  to  downstream  valve  closure 
was  investigated  for  the  KC-135  and  a  laboratory  test  rig  and 
compared  to  available  experimental  data  and  previous  research. 
Parametric  studies  were  performed  of  system  variables:  valve  closure 
time,  valve  area  versus  time  (closure  curves)^  and  variations  in 
surge  attenuation  components  (surge  boot,  accumulators),  A  typical 
closure  curve  for  quick  disconnects  was  approximated  and  verified. 
The  simulation  results  compared  favorably  with  experimental  data 
and  previous  work.  Transient  pressure  was  found  to  be  sensitive  to 
accumulator  or  surge  boot  precharge  pressure  while  changes  in 
accumulator  volume  had  little  effect.  Variations  in  accumulator 
entry  line  diameter  or  length  affected  transient  pressure  and 
settling  time.  Increased  valve  snubbing  reduced  maximum  transient 
pressures. 


DYNAMIC  CKARACTEHISTIC3  0?  AERIAL 
REFJELING  SYSTEMS 

I.  Introduction 

Problem  Rackfcround 

An  aerial  refueling  system  is  an  auxiliary  pumping  system  installed 
in  a  tanker  aircraft  which  is  used  to  transfer  fuel  to  a  receiver 
aircraft.  During  inflight  refueling  operations  fuel  is  pumped  from 
the  tanker  to  the  receiver  through  either  an  extensible,  rigid 
refueling  boom  or  a  flexible  hose.  Only  the  rigid  boom  method  will  be 
considered  in  this  study. 

Fuel  spillage  during  planned  or  inadvertent  separation  of  the 
tanker  and  receiver  is  prevented  by  use  of  a  spring  loaded  shut-off 
valve  at  the  boom  delivery  nozzle,  A  spring  loaded  valve  is  also 
provided  at  the  receiver  fuel  receptacle. 

The  rapid  closing  of  the  nozzle  shut-off  \’alve  during  separation 
can  cause  transient  disturbances  to  propagate  through  the  tanker 
refueling  system.  Depending  on  the  time  interval  of  the  transient, 
pressure  waves  traveling  at  speeds  close  to  that  of  sound  waves, 
about  3000  ft/sec,  can  easily  generate  energy  enough  to  damage  the 
system.  This  is  the  so-called  "waterhammer  effect".  Refueling 
system  performance  enhancements  such  as  higher  fuel  flow  rates  and 
faster  closing  shut-off  valves  for  less  fuel  spillage,  can  increase 
the  potential  for  damage  due  to  transient  pressure  surges.  Therefore, 
pressure  surge  effects  are  of  primary  interest  in  the  design  of  new 
systems  and  the  upgrading  of  existing  configurations.  For  the 


syster.  stru'.ture  to  be  maintained  within  aircraft  v/eight  standards 
and  still  withstand  "vrorst  case"  pressure  transients,  methods  must 
’oe  employed  to  attenuate  excessive  pressure  surges, 

'.'.’hen  designing  surge  attenuation  components,  reliable  design 
data  is  needed  to  make  a  reasonable  estim.ate  of  pressure  surge 
-lagnitude.  This  design  data  is  commonly  obtained  from  either 
analytical  studies,  experiment,  or  both.  General  fluid  transient 
phenomena  have  been  researched  and  documented  rather  extensively, 
’■.'ylie  and  Streeter  (Ref  l)  provide  a  comprehensive  study  of  fluid 
transients  and  Goodson  and  Leonard  (Ref  2)  summarize  the  development 
of  models  for  fluid  line  applications.  Streeter  and  Lai  (Ref  3) 
used  numerical  methods  and  the  digital  computer  to  solve  waterhammer 
problems;  their  solutions  were  verified  by  experiment  with  good 
results.  The  McDonnell  Douglas  Corporation  (Ref  4,5,6)  used 
computer  techniques  to  simulate  complex  aircraft  and  spacecraft 
hydraulic  systems.  Numerous  experiments  were  performed  to  verify 
various  test  cases. 

However,  for  aircraft  refueling  systems,  there  is  a  limited 
amount  of  well  documented  transient  response  experimental  data. 

Some  of  the  available  experimental  results  have  been  used  in 
extending  the  applications  of  existing  hydraulic  analysis  programs 
to  refueling  systems.  Parks  and  Franks  (Ref  ?)  investigated  modeling 
of  the  KC-135  tanker  system.  Their  study  focused  on  the  transient 
conditions  during  and  after  a  0,05  sec  closure  of  the  nozzle  shut-off 
valve.  Kinzig  et  al,  (Ref  8)  studied  the  com.patibility  of  the  KC-135 
and  KC-10  tankers  with  a  simulated  receiver  system.  Shut-off  valve 
closing  times  of  approximately  C,22  and  1,0  sec  were  used  to  simulate 


the  receiver  shut-off  valve  closures  during  refuelinsr.  The  objective  o: 


the  study  was  to  predict  transient  pressures  in  the  receiver  aircraft. 

Both  cf  the  refueling  system  simulations  discussed  above  agreed 
favorably  with  experimental  results  for  the  conditions  simulated. 
However,  there  is  an  element  of  uncertainty  critical  to  all  refueling 
system  simulations  that  must  be  resolved  to  achieve  accurate  solutions. 
This  element  is  the  shut-off  valve  closure  characteristics,  i.e.,  valve 
area  versus  time.  Kinzig  (Ref  8:10)  noted  that  small  changes  in 
closure  curve  shape  could  result  in  large  differences  in  pressure 
surge  magnitude.  A  McDonnell  Douglas  Corporation  (Ref  9*6-?) 
laboratory  test  indicated  nozzle  valve  closure  times  area  function 
of  separation  rate.  Thus,  for  varying  separation  rates  of  tanker 
and  receiver  the  valve  closure  characteristics  are  likely  to  change. 

As  a  consequence,  changes  in  separation  rate  will  influence  pressure 
surge  magnitude.  Therefore  an  accurate  approximation  of  value 
closure  characteristics  is  a  critical  input  for  refueling  system 
simulations. 

Scope  and  Objectives 

The  purpose  of  this  thesis  is  to  in\’estigate ,  by  computer 
simulation,  the  effects  of  variations  in  valve  closure  characteristics 
and  surge  attenuation  components  on  the  d;>Tiamic  response  of  an 
aircraft  refueling  system.  HYTRAK,  a  licDonnell  Douglas  Corporation 
(Ref  4)  developed  fluid  transient  analysis  computer  program  will  be 
used.  Two  distinct  system  configurations,  the  KC-135  and  a 
preliminary  KC-10  laboratory  test  set  up,  hereafter  referred  to  as 
the  Laboratory  Test,  will  be  modeled  and  simulated  with  the  program. 


The  specific-  objectives  of  the  study  are  to: 

a)  Establish  an  approx inatior.  of  the  actual  valve  closure 

4 

characteristics  for  a  typical  refueling  nozzle, 

b)  Verify  the  general  utility  of  HYTRAM  by  coTiparison  vfith  results 
of  Streeter  and  Lai, 

c)  Co.T.pare  the  si;nulatsd  results  of  this  study  with  experimental 
results  of  the  Laboratory  Test, 

d)  Determine,  for  the  Laboratory  Test  model,  the  influence  of 
variations  in  closure  curve  shape,  i,e,,  valve  area  versus  tine,  and 
modifications  to  surge  attenuation  components, 

e)  Determine,  for  the  KC-I35  model,  the  influence  of  both 
closure  curve  shape  variation  and  closure  time, 

f)  Examine  the  effect  of  varying  the  KC-135  surge  boot  precharge 


pressure. 


General 


Tanker  aircraft  are  equipped  to  transfer  fuel  to  receiver 
aircraft  by  means  of  an  extensible,  rigid  boom  or  a  flexible  hose 
arrangement  known  as  the  probe  and  drogue  system.  This  report 
considers  only  the  rigid  boom  method.  The  refueling  boom  is 
basically  a  telescoping  tube  attached  to  the  underside  of  the  fuselage 
that  is  maneuvered  into  contact  with  a  receiver  aircraft.  A  nozzle  at 
the  end  of  the  boom  fits  into  a  refueling  receptacle  on  the  receiver 
aircraft  to  provide  a  continuous  flow  path  between  the  two  aircraft, 

A  spring-loaded  poppet  valve  forms  the  seal  in  the  end  of  the  nozzle 
when  the  nozzle  is  not  in  contact  with  the  receiver  receptacle.  The 
poppet  valve  is  automatically  opened  during  insertion  of  the  nozzle 
into  the  receptacle  by  an  actuator  in  the  receptacle. 

The  primary  surge  attenuation  component  is  the  surge  boot,  A 
surge  boot  is  an  expandable,  gas  charged  envelope  that  surrounds  a 
perforated  section  of  the  boom,  norriallv  the  aft  portion.  The  basic 
parts  of  the  surge  boot  are  shown  in  Fig  1;  the  inner  surface  of  the 
envelope  seals  the  perforations  in  the  boom  until  the  fluid  pressure 
in  the  boor,  exceeds  the  precharge  or.  the  gas.  As  the  pre charge  is 
exceeded,  Aiel  under  pressure  forces  the  envelops  in.ner  surface  away 
from  the  cuter  surface  of  the  boom.  Therefore,  fuel  outflow  through 
the  perforations  and  into  the  cavity  compresses  the  gas  in  the  surge 
boot  and  provides  energy  dissipation  during  pressure  surges. 

During  a  normal  disconnect  of  tanker  and  receiver  a  signal  will 
be  automatically  sent  through  the  system  that  shuts  off  the  refueling 


Fig  1,  Surge  Boot  Details 


pumps,  opens  a  bypass  valve,  and  initiates  boom  retraction.  The 
opening  of  the  bypass  valve  allows  fuel  to  flow  back  to  the  tanker 
fuel  tank  so  there  is  no  sudden  flow  stoppage  during  disconnect  and 
retraction.  Any  surge  produced  under  these  conditions  will  be  mild  and 
easily  absorbed  by  the  surge  boot,  Honrial  disconnects  pose  little  or 
no  threat  to  the  system. 

The  abnormal  disconnect  conditions  that  cause  concern  are  the 
breaicaway  situations  due  to  turbulence,  emergency  actions,  or  a 
failure  of  the  automatic  process.  In  these  cases  it  is  likely  that 
there  will  be  a  flowing  disconnect  whereupon  the  pumps  continue  to  run 
against  the  closed  shut-off  valve  in  the  nozzle  with  no  retraction  of 
the  boom  or  bypass  of  the  flow. 

Failure  of  the  pumps  to  shut-off  increases  the  pressure  in  the 
line  to  deadhead  pump  pressure  (in  some  cases  90  psia  or  more)  in 
ad.dition  to  any  surge  pressure  induced  by  the  nozzle  shut-off  valve 


closure  (Ref  11:6,16).  Since  no  bypass  of  fuel  occurs  there  is  a 
stoppage  of  the  full  flow  through  the  boor:  at  the  tine  of  disconnect, 
■'■'igh  separation  rates  of  the  tanker  and  receiver,  up  to  the  nomal 
rn.  iinuir.  of  6,5  ft/sec  (Ref  9*7) »  are  usually  associated  with  abnormal 
disconnects. 


KC-136  System  (Ref  11,12) 

A  simplified  diagram  of  the  KC-135  refueling  system  is  shown  in 
Fig  2,  The  KC-I35  refueling  system  is  a  four  pump  (constant  rpm, 
each  rated  at  JOO  gum  0  80  psig)  system  with  two  refueling  pumps  in 
the  forward  body  tank  and  two  in  the  aft  body  tank.  Reverse  flow  is 
prevented  to  each  pump  by  individual  check  valves.  Fuel  flows  from 
all  four  pumps  are  united  in  a  main  trunk  line;  fuel  then  flows 
through  a  shut-off  valve,  pressure  regulator,  and  venturi.  Finally 
the  fuel  flows  through  the  boom  which  contains  a  surge  boot  with  a 
nominal  precharge  of  50  Psig*  The  pressure  regulator  is  set  to 
maintain  the  steady  state  pressure  at  the  nozzle  to  50  “  5  Psig. 

A  secondary  method  of  pressure  surge  protection  is  provided  by  two 
level  control  valves  located  in  the  lines  leading  from  the  refueling 
manifold  to  the  forward  and  aft  body  tanks.  The  valves'  cracking 
pressure  of  S5-115  psi  allows  fuel  to  flow  to  the  body  tanks  during 
pressure  surges  (Ref  10:3). 

Laboratorrr  Test 

The  McDonnell  Douglas  Corporation  performed  an  experimental 


investigation  (Ref  9)  to  study  surge  pressures  developed  by  the 
closure  of  shut-off  valves  during  tanker/rcceiver  disconnects.  These 
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Fig  3»  Lauoratory  Test  Configuration 


Ill,  Con-.puter  ?ro:/ra.T.  Ziescription 

General 

HYTRAII,  the  hydraulic  transient  analysis  computer  program  used  in 
t'lis  study,  was  developed  by  the  KcDonnell  Douglas  Corporation  under 
an  Aero  Propulsion  Laboratory  crntract  (Ref  4),  The  contract  included 
not  only  program  development  but  also  extensive  \'alidation  and 
documentation.  The  program  uses  a  building  block  approach  to 
facilitate  simulation  of  complete  systems.  After  the  configuration  is 
established  according  to  HYTRAN  format,  the  program  calculates  the 
values  of  flows,  pressures,  and  component  variables  throughout  the 
system. 

The  basic  theory  of  HYTRAK  is  that  a  change  in  flow  (momentum) 
results  in  a  change  in  pressure,  Plow  is  brought  to  rest  at  the 
expense  of  an  increase  in  pressure,  3y  coupling  the  continuity  and 
momentum  equations  with  component  equations,  pressure  and  flow  as  a 
function  of  time  and  line  distance  can  be  determined.  The 
mathematical  solution  technique  used  by  the  program  is  the  method  of 
characteristics.  The  major  assumptions  upon  which  the  program  is 
based  are; 

a)  Fluid  temperature  is  held  constant  during  the  entire  run, 

b)  Flow  is  one  dimensional,  that  is,  the  fluid  properties  are 
constant  across  any  transverse  cross  section  of  pipe, 

c)  Pipes  have  circular  cross  sections  and  friction  factors 
based  on  smooth,  drawn  tubing, 

d)  Stresses  In  pipes  are  always  below  the  elastic  limit. 


g8c^ie^r^•  is  such  that  the  "thi::  case  is  •.•aiiu. 


f]*  Pipe  ar.i  lieuid  are  perfectly  ei:-e*.ic  (^.11  er.erry 
dissipation  is  due  to  shearing  stresses  a*,  tr.e  •.■alls]' , 

Each  component  subroutine  has  its  ovm  set  of  as.svnptior.s  since  each 
model  is  broken  dov/n  into  its  most  basic  eouaticr.s  :'f  .".otion  and  flo 

The  controlling  input  to  the  system  is  ncr^'aHy  a  -v-al’/e  notion, 
i-jhich  causes  a  disturbaiiCe  to  propagate  thr:‘-gh  nhe  oyster;  at  or  nea 
the  speed  of  sound.  The  components  of  "ohs  system  each  ha'/e  a 
response  to  the  pressure  and  flow  changes. 

HYTP^N  is  composed  of  four  basic  parts:  input,  steady  state 
calculations,  transient  calculations,  and  output.  A  simplified 
flowchart  is  shown  in  Fig  4, 

Included  in  the  input  data  are  line  and  component  parameters, 
the  system  configuration,  i.e,,  how  the  lines  and  components  are 
connected,  initial  conditions,  fluid  properties,  and  various 
directive  parameters  such  as  tine  steps,  run  time,  plotting 
intervals,  etc. 

The  steady  state  section  of  the  program  balances  the  pressures 
and  flows  in  the  system  and  calculates  the  initial  values  for  all 
the  system  \'3,ri?bleE,  Cnee  the  initial  %’alues  are  established  at 
zero  time,  the  program,  starts  by  calculating,  fer  a  small  change  in 
tine  (at),  new  flows,  pressures,  and  values  for  the  component 


variables  s'uch  as  valve  position,  accum'. 


Luid  volu.-.es,  etc. 


The  program  continues  forward  in  AT  time  ir.ter.'als,  first 
calculating  the  line  and  then  the  component  •.■•ariables  until  the 
desired  simulation  time  period  is  covered. 

The  output  is  time  histories  of  seisemei  system  variables  which 
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Fig  4,  Simplified  HYia^N  Flowchart 


have  been  a _ -.ur’xi  by  the  controlling  input,  A 

V,'C'l*^d  06  c.  DZ'SSSU^T^  VCiTSMC  “tinS  CUlT'V'r  f'oiT  C6l'0X*£ 

a  valve  closure . 


A  detailed  description  of  pa-cgrar’  str'ucture 
ana  conponent  node Is  is  provided  in  Ref  4,  The 
approxir.ately  170,000  core  nenory  and  a  r’un  tine 


typical  output  plot 
during,  ana  after 


netheds  of  solution 
TOgra:.  requires 
betv.-een  80  and  300  sec 


Fuel  Systen  Usage 

Although  HYTPAI'  was  develoi^ed  for  use  with  hydraulic  syster.s,the 
basic  assunptions  and  solution  techniques  are  corjnon  to  most  fluid 
systems.  The  main  alterations  required  for  use  with  fuel  systems 
was  to  change  the  fluid  properties  fron  hydraulic  fluid  to  fuel. 


System  modeling 

The  first  step  in  developing  a  HYTRAN  model  is  to  construct  the 


system  scnematic  in  HYTRAN  formiat.  This  involves  the  assignment  of 


line,  component,  node,  and  leg  numbers  throughout  the  system  as  well 
as  determining  which  component  models  will  be  selected  to  simulate 


each  individual  component.  References  13  and  14  explain  in  detail  how 
to  translate  a  system  into  HYTRtN  format. 

The  specific  HYT?u\H  component  models  used  for  the  simulation 
miodels  discussed  in  Chapter  IV  and  V  are  discussed  briefly  in  the 
following  paragraphs.  Detailed  descriptions  of  each  component  model 


can  be  found  in  Ref  4  and  15, 

Fuel  tanks  and  flow  sources  were  simulated  with  HYTRAN  TYFZ  6l 
constant  pressure  reservoirs.  The  reservoir  is  assumed  to  have  an 


infinitely  large  gas  volume  so  that  pressure  remains  unchanged 


Lines  were  joined  at  branches  and 


'■zior.c  wizn 


branches.  The  branches  do  not  incoroor?-t£  any  losses  associated 
chanf^es  in  diameter  and  flow  direction  or  divisions,  i.e,,  they  are 


frsctionless  branches. 


Fixed  orifices  were  modeled  using 


TiPE 


oostricncrs , 


The  model  asaur.es  flov:  in  either  direction 
is  the  sane  for  flow  in  each  direction  and 


,  The  discharge  coe 
constant  over  the  entire 


regime , 

The  FYTRAN  TYr2  31  check  v'alve  model  was  used  to  sir.ulate  chec 
and  relief  valves.  The  check  valve  is  r.oieled  as  an  orifice  of 
variable  area.  The  area  at  any  given  time  is  determined  by  the 
position  of  the  force-balanced  poppet.  Reverse  flow  can  take  place 
until  the  valve  closes. 

Control  valves  were  simulated  using  FYTRAN  TYP2  21  tvro-way 
Nralves,  The  val\'e  is  modeled  as  an  orifice  of  -/ariablc-  area  and  a 
constant  discharge  coefficient  over  the  entire  flew  revnge ,  A  time 
history  of  valve  closure  is  input  in  the  form  of  an  effective  valve 
area  versus  time  table.  The  effective  valve  is  defined  as  the 
valve  area  at  a  particular  ti.mc-  multiplied  by  the  discharge 
coefficient,  In  this  study',  valve  area  versus  time  plots  are 
presented  in  the  r.or.d'’~.ensicnal  form  of  val'/e  position  versus 
percent  of  closure  cine.  The  valve  pc.eitio.'. ,  given  in  percent  open 
is  the  valve  ar«a  with  respect  to  the  .-..Id  flow  valve  area,  Li.':ev:i 
the  percent  of  closure  time  values  are  -.he  time  '/alues  for  each 

r 

valve  position  with  respect  to  the  total  closure  time,  Therefore, 
the  effective  area  w.' th  respect  tc  ti.me  for  prograr.  input  is  obteadn 
by  .multiplying  the  full  flow  area  by  valve  position  and  discharge 


coefficient 


-3  choun  celow 


A 


K(t)C^A 


VO 


(1) 


vjhere 


A  (t)  =  effecti'.’e  valve  area  '.rith  resoect  to  time 

K(t)  =  valve  position  with  respect  to  time 

C,  =  discharge  coefficient 

d  ^ 

A  =  full  flow  valve  area 

VO 


Accumulators  were  simulated 
Each  unit  is  modeled  as  a  simple 
with  the  frictional  and  inertial 
a  nonvarj’’ing  gas  constant  and  no 


with  KYTRAK  TYPE  ?!  accumulator  models, 
gas  charged  piston  type  accumulator 
effects  neglected.  The  model  assumes 
heat  transfer  between  the  gas  and  the 


accumulator  walls 


Laboratory  Test  Cor.ou'ter  iiodel 


The  Laboratory  Test  configuration  in  Fig  3  modeled  using  the 
HYTRAN  conponer.ts  discussed  in  Chapter  III.  A  schematic  of  the  test 
configuration  in  KYTRAF  input  format  is  shown  in  Fig  5* 

Input  data  for  all  lines  except  those  entering  the  accumulators 
vrere  obtained  from  the  Laboratory  Test  report  (Bef  9)  •  The  outside 
diameter  of  the  accumulator  entry  lines  was  not  available  and  was 
estimated  from  measurements  of  test  rig  photographs.  An  assumed 
inside  diameter  was  based  on  the  thickness  of  other  tubing  used  in 
the  experiment. 

Fluid  properties  were  altered  from  hydraulic  fluid  to  those  of 
JP-5  at  60  F,  The  properties  used  were: 

Viscosity  =  0,0362  X  10”^  in^sec 

Density  =  0,?20  x  10  ^  Ib^-sec^/in^ 

Bulk  Modulus  =  0.139  x  10"  Ib^/in 

Vapor  Pressure  =  0,5  psia 

Time  steps  ( AT)  of  0,0005»  0,0009-,  and  0,000333  sec  were  used  over  a 
simulation  interval  of  0,4  sec.  This  interval  was  chosen  to  match  the 
time  period  over  which  the  expierimental  data  were  recorded.  The 
various  tine  steps  were  used,  when  required,  to  insure  a  whole  number 
of  time  steps  over  all  segments  of  each  of  the  effective  valve  area 
versus  tine  input  tables. 

The  nozzle  was  modeled  using  a  TYPE  21  two-way  valve.  The  valve 
closure  was  input  as  an  effective  \'alve  area  versus  time  as  discussed 
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TYPE  71  Aceuaulator 


— >  TYPE  6l  Constant  Pxssaure 
Reservoir 

—  TYPE  11  Branch 


~  TYPE  21  Two-saj  valve  with 

nodes  at  each  connection  j, 

— •  Node  nusher  . 

—  Conponent  nusher  [ 

—  Line  nusher  ^ 


Fig  5«  HYTBAN  Diagras  of  Laboratory  Test 


in  Chapter  III 


.  This  method  sinplified  nodifications  and  the  usare 
cf  '.'arious  closure  our/cs. 

The  full  flow  ai«a  tines  discharge  coefficient  was  calculated 
from  a  flow  versus  pressure  difference  across  the  nozzle  cur'/e  for 
the  actual  nozzle  (Ref  16:44),  A  flow  rate  of  2016  in''/ sec  '^as  use 
in  the  Laboratory  Test.  This  flovirate  corresponds  to  A?  =  $,1  csi 
from  the  curve  for  the  nozzle.  Using  the  orifice  equation, 


2 

a  steady  state  value  of  =  5*35  in  was  computed. 

The  steady  state  flow  was  exhausted  into  a  TYPE  61  constant 
pressure  reservoir  to  simulate  the  fuel  return  system  used  in  the 
test.  The  constant  pressure  was  set  at  4-7,5  psia.  This  pressure 
was  computed  by  subtracting  the  AP  of  the  nozzle  from  the  measured 
steady  state  flowing  pressure  of  52.6  psia  just  upstream  cf  the 
nozzle  (Node  ll) .  This  back  pressure  was  consistent  with  the 
conditions  used  to  generate  the  flow  curve  for  the  nozzle, 

TYPE  71  gas  charged  accumulators  v;ere  used  to  simulate  the 
test  accumulators.  Volume  and  precharge  pressure  was  the  only 
definitive  data  available  for  the  accumulators;  therefore,  some 
assumotlcns  were  required  to  complete  the  models.  Sntr^/  lines  were 
scaled  from  photographs  to  be  6  in  long  with  an  outside  diameter  of 
2  in,  A  wall  thickness  of  0,125  in  wa-s  assumed.  The  volume  of  sac 
accumulator  was  given  as  5  cal  (1155  in  ) »  a  usable  volume  of 

•3 

1000  in-'  -was  used  in  the  simulation  to  account  for  the  unknown 


internal  cc:  ''i.-uration. 


AssuHiing  that  the  ~as  v;ould  be  cor.pressei 
to  a  r.axir.'jn  of  ten  percent  of  its  original  volume  resulted  in  a 
na-xinum  fuel  volume  for  each  accumulator  of  900  in^.  The  total 
fuel  volume  for  the  four  accumulators  totaled  3'^00  consistent 

with  the  actual  volume  of  35-5  (Ref  1?)  for  the  KC-10. 

!'o  data  other  than  photographs  was  available  for  the  Laboratory 
pumping  system.  Since  the  system  provided  a  constant  smeady  state 
flow  and  pressure,  a  TYPE  6l  constant  pressure  reservoir  was  used  to 
model  the  pump  as  a  constant  pressure  source.  The  pressure  of  the 
reservoir  was  adjusted  to  provide  the  test  condition  steady  state 
flow. 

A  complete  listing  of  input  parameters  for  the  Laboratory  Test 
model  is  given  in  Appendix  A, 


The  TT-I35  systen  was  modeled  usin? 
(Ref  10)  and  the  r.YZ3I\'.l  component  models 
The  use  of  a  similar  model  and  test  condi 
and  extension  of  the  orevious  results,  A 
ea-rlier  version  of  uYTRAT  which  required 
parameters,  the  basic  model  configuration 


,e  technique  of  ?a,rl;s 
3 sussed  in  Chapter  III. 
cns  provided  a  cross-check 
,hou:.'h  Farhs  used  an 
0  or. nutation  of  some  input 
■as  unchanged.  A  schematic 


of  t: 


he  system  in  "YTRAM  format  is  shown  in  Fig'sres  6  and  ?, 
Fluid  properties  were  altered  from  hydraulic  fluid  to  JP-4 


at  60  F,  The  properties  used  wares 


Viscosity 
Density 
Bulk  modulus 
Vapor  Pressure 


0,136  X  10  ^  in^/sec 
0,77  X  10  ^  Ib^-sec^/in^ 
0.16  X  10"^  IbVin^ 

2.0  psia 


Time  steps  (AT)  of  0,0005,  0.0004,  and  0,000'}'}’}  sec  were  used 
over  a  simulation  interval  of  1,0  sec.  This  interval  was 
selected  to  match  Parks'  simulation  time  innsr/al  (Ref  10j23). 

The  various  time  steps  were  used,  when  required,  to  insure  a  whole 
number  of  time  steps  over  all  see:nents  cf  the  •'/alve  closure  input 
tables. 


Parks  simulated  the  KC-135  system  with  three  pumps  on  and  one 
pump  off  (?;ef  10:23);  this  study  used  the  same  configuration  to 
facilitate  ccr.pariscns  with  Parks.  Zach  nur.p  was  mcieled  as  a  TYPE 
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Fig  6.  HYTRAN  Schenatlc  of  KC-133  Refueling  System 
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Pig  7«  HYTRAN  Schenatlc  of  KC-135  Surge  Boot 


7^ 


T-TT^rt 


were  sinula^_-d  as  TV?2  31  checlr/alves. 

The  line  valve  and  re-^ulator  were  nodeled  as  TY??  ^1  restrictors; 

the  level  control  valves  as'TYPT  3I  checltvalves  and  the  fuel  tanks  as 

TYPE  61  reservoirs,  respectively.  The  surge  boot  was  nodeled  as 

twelve  discrete  accumulators;  the  total  volume  of  the  twelve  units  was 
3 

1300  in-"^.  Each  was  modeled  as  an  individual  TYPE  71  accumulator  with 
a  minimum  gas  volume  equal  to  ten  percent  of  the  initial  gas  volume, 

A  TYPE  21  control  valve  was  used  to  simulate  the  boon  nozzle. 

The  valve  closure  was  input  as  an  effective  -.'alve  area  versus  time  as 
discussed  in  Chapter  III. 

All  dimensions,  where  applicable,  were  obtained  from  the  Parks 
model  (Ref  10 1 70-77). 

The  steady  state  flow  was  exhausted  into  a  TYPE  6l  reservoir  to 
simulate  the  fuel  return  system.  The  reservoir  pressure  was  set  to 
43  psig  (Ref  10:36). 

A  complete  listing  of  input  parameters  car.  be  found  in 
Appendix  B, 
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VT 
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Disci!3sic,r.  cf  Valve  Closure 


General  (Ref  l';69) 

As  discussed  in  Chapter  1,  a  r.czzle  valve  closure  causes  a 
pressure  surj^e  to  propagate  through  a  refueling  systen.  The  direc 
relationship  of  pressure  surge  nagnituds  tc  valve  closure  can  he 
illustrated  by  considering  the  deceleraoion  of  an  inconpressihle 
fluid  flowing  frictionlessly  in  a  rigid  pipe  of  unifom  area  A  v;it: 
a  velocity  V.  The  pipe  has  a  length  L,  an  inlet  pressure  ,  and 


a  pressure  P2  at  L.  Assume,  at  length  L,  there  is  a  \'al'/e  'which  c 


reduce  the  velocity  at  L  to  V  -  A7,  The  mass  rate  of  flow  for  a 
pressure  -wave  traveling  at  sonic  velocity  c,  is  r.  =  yAc,  From  th 
impulse-momentum  equation,  for  this  application 


(  ^Ac)  (7  -  AV  -  7)  =  P-A  -  F^A 

or  the  increase  in  pressure  is  given  by 

AP  =  -  jc  AV  ( 

thus  a  change  in  velocity  is  directly  proportional  to  a  change  in 
pressure.  The  tine  for  a  press'ure  'wave  to  travel  the  length  of 
pioe  L  and  return  is  t  =  2L/c,  If  the  time  of  closure  T  is  less 
than  or  equal  to  t,  the  approximate  pressure  rise  is  given  by 

AP  -  -2:3 

‘C 


thus  a  deer’  .36  In  closure  time  will  cause  an  increase  in  pressure 
rise , 

Since  the  reduction  in' flow  velocity  is  determined  by  the 
valve  closure ,  it  is  apparent  that  the  closure  curve  is  a  critical 
input  to  the  simulation  of  aircraft  refueling  system  response 
subsequent  to  nozzle  valve  closures. 

Valve  Closures  from  Previous  Research 

Parks  (Bef  10:19-20)  used  a  val’/e  closure  relationship, 
hereafter  referred  to  as  curve  "A”,  described  by  Eq  (l)  where 


The  closure  curve  corresponding  to  this  relationship  is  shown  in 
Fig  8  and  is  characterized  by  a  rapid  reduction  of  the  valve 
opening  initially,  followed  by  an  almost  asymtotic  reduction  to  full 
closure c  This  slower  area  reduction  during  the  final  stage  of  the 
closure  curve  is  referred  to  as  "valve  snubbing".  Parks  justified 
his  usage  of  the  curve  on  the  basis  of  a  compearison  to  experimenta’ 
data. 

Typical  curves  used  by  Kinzig  (Ref  8:25)  in  his  study  of  receive 
pressure  surges  are  shown  in  Fig  9» 

For  use  in  this  study,  the  curves  were  reduced  to  an  analytical 
expression  of  the  following  form: 


Fig  8,  Closure  Curve  *A' 


i:(t)  = 

fl  -  (t  /  T^) 

J  [a  -  b(l  - 

^  (n) 

(t  / 1^)3 

ere  for 

curve  "5” 

curve  "C" 

a  =  1.45 

a  =  2.60 

b  =  l.OC 

b  =  2.00 

n  =  1.30 

n  =  1.25 

th  curves  represent  a  relatively  gentle  closure  characteristic 


Fig  9«  Closure  Curves  *B*  and  *C 
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(Rof  and  drawings  of  the  valve  nechanisr.  were  used  to  .olot 

actual  closure  curve  fra'  a  rerxielinr  nonule.  The  curve,  in 
nondr! ueneionr.l  for',  is  shown  in  Pir  1C,  Pror'i  fully  open,  th 
closes  .it  an  increasi.nr  rate  as  it  accelerates  to  a  constant  closure 
S'ate,  The  closure  rate  decreases  during  the  final  stage  of  area 
reduction,  approximately  the  last  ten  percent.  This  snubbed  closure 
is  expected  since  the  \'alve  poppet  incorporates  a  damper  that  is 
designed  to  slow  the  closure  during  the  last  stage  of  poppet  travel. 


Nozzle  /  Recentacle  Combination 

VJhile  the  previous  curves,  Figs  8  and  5,  may  be  tj'pical  for  normal 
disconnects,  abnormal  disconnects  at  high  separation  rates  are  likely 
to  result  in  more  severe  closures.  Before  establishing  the  effect  of 
high  separation  rates  or.  the  \’alve  closure  cur'/e,  the  interaction  of 
the  combined  nozzle  and  receptacle  config-uiation  we^s  examined, 

"hown  in  Fig  11  are  (a)  the  receptacle  sleeve  v-alve,  (b)  the 
nozzle  \'alve,  (c)  the  nozzle  and  receptacle  in  connected  position. 

The  spring-loaded  nozzle  poppet  valve  is  opened  .mechanice.lly  on 
contact  with  the  fixed  pedestal  in  the  receiver  aircraft  receptacle 
in  conjunction  with  the  cpeni.ng  of  the  sliding  sleeve  valve  in  the 
receptacle  by  the  nozzle  tip. 

During  quick,  disconnects,  the  damper  in  the  poppet  can  slow  the 
poppet  closure  such  that  the  poppet  and  pedestal  do  not  remain  in 
contact  during  separation  of  the  nozzle  and  receptacle,  Th.ereiore, 
nozzle  valve  closure  is  independent  of  separaiior.  rate.  Conversely, 


..r 


Fig  10 


,  Nozzle  Closure  Curve  Coiriputed  from  the 
Laboratory  Test  Neasurements 

since  the  receptacle  sleeve  valve  is  spring-loaded  and  undamped,  its 
closure  rate  is  determined  by  the  rate  of  separation.  A  typical  closure 
curve  for  the  receptacle  sleeve  ’/alve,  hereafter  referred  to  as  cur-v'e 
"D" ,  is  shcun  in  Fig  12.  The  closure  accelerates  to  a  constant  rate 
with  no  appreciable  snubbing. 

To  reach  the  fully  closed  position,  the  required  travel  for  the 
receptacle  sleeve  \'alve  is  1.0  in  (Ref  9*7)  and  for  the  nozzle  valve 
poppet  a  minimum  of  l.lp  in  is  required  to  reach  the  closed  position 
(Ref  I6t50).  Therefore,  when  the  nozzle  and  receptacle  have  separated 
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(a)  Receptacle  Sleeve  Valve 


1.15  in 


Spr?.n5-loaded 
Poppet  - 


^  \ 
\  S 


Poppet  Valve  in 
closed  Dosltion 


Poppet  Valve  in 
open  position 


(b)  Nozzle  Valve 


Fig  11.  Nozzle  /  Receptacle  Details 


Fig  13.  Nozzle  and  Receptacle  at  the  One  Inch 
Separation  Point 

system  even  though  the  nozzle  valve  is  still  open.  This  supposition 
is  supported  by  pressure  measurements  obtained  during  the  Laboratory 
Test  (Hei  9!l0)  which  indicate  peak  pressure  occurs  coincident  with 
the  1.0  in  separation  point  for  separation  rates  of  f-lO  ft/sec. 


Actual  Valve  Closure  Anproximation 

V/ith  a  relationship  between  separation  rate  and  valve  closure 
established,  a  reasonable  approximation  tc  the  r.inir.um  effective 
closure  time  can  be  established.  The  r.axir.ur.  disconnect  speei  is 
6.5  ft/sec  (Ref  9«7).  Assuming  instantaneous  acceleration  tc  this 


rate ,  the  time  required  for  the  nozole  and  the  receptacle  to  separate 
one  inch  is  0,013  sec.  The  separation  rates  measured  during  the 
Laboratory  Test  (Ref  9!l0)  indicate  a  r.inimur.  time  for  one  inch 
separation  of  approximately  0,02  sec  at  a  separation  rate  of  6  ft/sec 
Based  on  the  previous  discussion  in  this  chapter  the  closure 
curve  for  quick  disconnects  is  theorized  to  be  basically  the 
receptacle  closure  curve  as  shorn  in  Pig  12  (curve  D),  Typically  the 
closure  rate  would  increase  to  a  constant  rate  vfith  no  snubbing.  The 
overall  effect  of  the  combined  nozzle/receptacle  closure  is  to 
eliminate  the  snubbing  portion  of  the  closure  action  of  the 
receptacle , 


Vllt  fiesults  and  Di:^cucsiori 


IIYTR4N  Verification 

A  basic  Katerharriiner  problen  with  known  results  vras  solved  using 
HYTRAN,  Streeter  and  Lai  (Hef  3)  obtained  ar.  analytical  solution  for 
the  configuration  shown  in  Fig  14  and  verified  it  by  experiment. 


After  ertaclishing  a  steady  state  flowrane 
closed  according  to  Sq  (l)  where 


the  line ,  the  valve  vras 


K(t)  = 


C7) 


The  HYTHAN  format  configuration  is  shown  in  Fig  15#  A  closure  time 
of  0.09  sec  and  a  steady  state  flowrate  of  5  in'' /sec  was  used.  The 
working  fluid  was  water  at  60  F. 

Comparison  of  KYTRAN  and  Streeter  results  (Fig  I6) ,  indicate 
excellent  agreement  in  phasing,  curve  shape,  and  magnitude  of  the 
first  peak.  Subsequent  peaks  of  the  HYTHAR  solution  exhibit  more 
attenuation  than  the  Streeter  solution  due  to  the  inclusion  of 
djmamic  friction  effects  in  HYTR/iN  that  Streeter  neglected. 


Laborator:/'  Test  Simulation 

Comparison  with  Experimental  Results,  The  steady  state 
conditions  prior  to  the  sepaiution  of  nozzle  and  receptacle,  i.e,, 
valve  closure,  vrere  as  described  in  Chapter  IV,  The  accumulator 
precharge  was  set  at  9C  psia  (Ref  9*10#  2L),  The  '.^alve  was  closed 


Fig  l6,  Compa^'ison  of  Transient  Response  Due  to  Valve  Closure; 
HYTRAN  versus  Streeter 


according  to  Hq  (l)  with  K(t)  from  curve  D  ci  Chapter  VI,  Pig  12. 
Closure  tiac  was  O.O3  oec.  An  oscillograph  araee  of  -.'.a;  •..-•.p&ri-.ental 
response  for  the  conditions  of  the  sinulaticn  was  obtained  frcn  the 
test  report  (Ref  9S24), 

Preliminary  results  exhibited  reflected  peaks  net  shewn  in  the 
exoerir.ental  data  that  indicated  the  laboratory  pumpiny  :y.-ster  did 
not  reflect  pressure  waves  precisely  as  a  constant  pressure  reservoir 
As  a  result,  the  laboratory  pumping  system  was  modeled  as  a  constant 
pressure  reservoir  connected  to  a  long  line,  the  length  of  which 
prevented  the  return  of  the  reflected  wave  during  the  time  interval 
of  interest  (0,0  to  0,4  sec).  The  line  from  source  to  test  section 
(LI)  was  lengt’nened  to  accomplish  this  modification  of  the  model. 

The  pressure  of  the  constant  pressure  source  was  adjusted  to  provide 
the  desired  steady  state  flow  conditions. 

Comparison  of  experimental  and  simula+ei  transient  response  to 
the  valve  closure  is  shown  in  Fig  17,  The  first  peak  in  both  cases 
corresponds  to  the  complete  closure  of  the  \-alve.  Excellent- agreemen 
is  seen  in  comparison  of  the  five  percent  settling  tines.  The 
settling  time  is  the  time  required  for  the  response  to  decrease  to  a 
specified  percentage  of  its  final  ^'alue,  The  ir.itial  and  final 
values  are  virtually  the  same  for  both. 

Comparison  of  peak  raagritudes  and  phasing  are  less  favorable. 

The  computer  simulation  underpredicts  the  magnitude  of  the  first  peak 
and  cverpredicts  the  subsequent  peaks.  There  are  twice  as  many  peaks 
in  the  experimental  results  as  in  the  simulation;  a  result,  perhaps. 


of  the  accumulator  model  limitations 
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Fij?  17.  Comparision  of  Experimental  ar.i  Sir.ulated  Transient 
Response  at  Nozzle  to  Valve  Closure  for  the 
Laboratory  Test  ilodel 


Accumulator  Volume  and  Precharfte,  Four  accumulator  volume 
confirpirations  •■:1th  four  different  precharr^e  sot'-lngs  per  configuration 
were  investigated.  All  simulation  conditions  v;ith  the  exception  of 
accumulator  parameters  were  outlined  in  the  previous  section.  Shown 
in  Figs  lS-21  are  the  maximum  transient  pressures  -/ersus  precharge 
pressures  for  each  volume  configuration.  Refer  to  Fir  3  for  the 
relati’/e  locations  of  the  boom  inlet  and  nozzle. 


—  Peak  Nozzle 


18.  Effect  of  Accumulator  Pxecharge  on  Maximum  Transient 
Pressure  for  Accumulator  Volume  of  1000  cubic  Inches 
for  the  Laboratory  Test  Model 


—  Peak  @  Nozzle 


Pig  19.  Effect  of  Accumulator  Precharge  on  Maximum  Transient 
Pressure  for  Accumulator  Volume  of  660  cubic  Inches 
for  the  Laboratory  Test  Model 


Fig  20.  Effect  of  Accumulator  Prechaxge  on  Maximum  Transient 
Pressure  for  Accumulator  Volume  of  $00  cubic  Inches 
for  the  Laboratory  Test  Model 


®  Nozzle 


Fig  21,  Effect  of  Accumulator  Precharge  on  Maximum  Tranolent 
Prescure  for  Accumulator  Volume  of  250  cubic  Inches 
for  the  Tjaboratory  Test  Model 


Volume  changes  had  a  negligible  effect  cn  t:.e  r.axl.'iUni  transient 
pressure  at  either  the  nossle  or  boo;,  inlet,  A  small  change  in  the 
second  peak  pressures  at  low  precharge  as  the  volumes  are  decreased  is 
due  to  the  accumule-tors  reaching  their  fluid  capacity  limits.  The 
dominant  influence  was  precharge  pressure,  increasing  pr^charge 
pressure  resulted  in  higher  peak  pressures  at  both  nozzle  and  boom 
inlet,  hete  that  the  first  peak  is  unaffecaed  by  changes  in 
either  volume  or  precharge.  This  is  due  to  the  accumulator  time 
response  lag.  The  initial  impulse  is  not  attenuated  as  readily  as 
surges  that  build  up  over  a  longer  period  of  time  since  a  finite 
period  of  time  is  required  for  flow  into  the  accumulators  to  begin. 
Accumulator  Entry  Line  Diameter,  Simulations  with  accumulator 
entry  line  diameters  from  1,25  to  2,75  in  ^nin  to  observe  the 


general  effects  of  diameter  variation.  Also,  since  the  actual 
accumulator  entry  line  diameters  for  the  Laboratory  Test  were  not  known 
a  verification  of  the  choice  of  1,75  in  was  desired.  All  simulation 
conditions  except  for  entry  line  diameter  were  as  outlined  in  the 
previous  section.  The  basis  for  comparison  was  the  maximum  transient 
pressure,  five  percent  settling  time,  a.nd  the  number  of  peaks  to 
settling  time. 

Changes  in  diameter  did  not  affect  the  .-.aximum  pressure  at  the 
boon  nozzle.  As  in  the  previous  section,  this  is  attributed  to  the 
accum.ulator  response  ^ime,  The  maximum  pressure  at  the  boom  inlet  is 
affected  because  during  the  time  required  for  the  pressure  wave  to 
reach  the  boom  inlet,  flow  into  the  accumulators  has  begun.  The 
increase  in  pressure  magnitude  at  the  boon  ir.let  for  a  decrease  in 
the  diameter  is  shown  in  Fig  22,  This  increase  is  expected  due  to  the 


for  the  Laboratory  lust  fiolol 


additioricil  resiatance  to  flow  in  snaller  alias,  thus  necaiing  sona  of 
the  capability  of  the  acou:nulator  to  function  as  onerry  ahnerplion 
component , 

Settling  time  and  number  of  peaks  to  settlinfg  time  were 
favorably  affected  by  decreasing  the  line  diameter,  ?ig  23  £-nd  2b, 
Settling  time-  was  reduced  by  2C  to  50  percent  for  the  range  of 
diameters.  The  nuinber  of  peaks  from  initiation  to  settling  tine  was 
reduced  ^0  to  ?C  percent.  In  this  case  the  increase  in  flow 
restriction  to  the  accumulator  resulted  in  increased  damping  and  less 
pulsation  added  to  the  system  due  to  the  release  of  energy  stored  in 
the  accumulators, 

Tv7o  runs  were  made  with  longer  than  standard  (14  in  longer) 
accumulator  entry  lines.  The  results  for  the  longer  lines  exhibited 
the  same  response  trends  as  decreasing  the  diameter,  see  Fig  23. 

Exam.ination  of  settling  tines  at  the  nozzle  supports  the 
choice  of  the  1.75  in  entry  line  for  the  Laboratory  Test  simulation. 
The  settling  time  of  approximately  0,02  sec  agrees  with  the  measured 
response  (Fig  1?) , 


Closure  Curve  Study 

Laboratory  Test  I.odel,  Simulations  were  run  using  each  of  six 
closure  curves,  shown  in  Fig  25,  to  examine  the  effects  of  curve 


shape.  Steady  state  conditions  were  as  described  in  Chapter  IV 
with  an  accumulator  precharge  of  90  psia.  Valve  closure  was 
according  to  Sq  (l)  vrith  K(t)  given  by  the  closure  curves.  Time  of 
closure  was  0,03  sec. 


iitrv  Line 


iiffoct  ot'  Variation  of  Accunulator  EntT7  Lino  Diameter  on  Trr.nr.iont  Hospo: 
i<'o'«nle  for  t’lO  Laboratory  Tcot  i.oiicl 


Accumulator  Entry  Line  Diameter  (in) 


Valve  Posi 


Fig  25*  Closure  Curves  for  Laboratory  Test  Model  Closure 
Curve  Study 
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TABLE  I 


;:.ifect  of  G’ar'v'e  Type  on  Transient  rrecsure  and 
Settling  Tine  for  the  Laboratory  Test  i'.oiel 


i-.axir.u:;.  .ransient  rre 


Gsara  (nsia) 


Curve 

Foszle 

Boom  Inlet 

Settling  Time  (sec; 

Linear 

3iC 

137 

.2^3 

D 

3^c 

210 

.246 

Z 

3^0 

190 

.232 

F 

340 

206 

.247 

G 

340 

220 

.247 

H 

260 

122 

.242 

Results  for  the  six  curves  are  shown  in  Table  I,  Curves  D,  E,  F, 
and  G,  produced  similar  results,  not  unexpected  since  the  general 
curve  shape  is  the  same  for  all  four  (Refer  to  Fig  25) #  Maximum 
transient  pressure  at  the  nozzle  and  boom  inlet  for  those  four  curves 
are  the  same  except  for  cui-ve  Z,  Curve  S  has  a  less  severe  closure 
vrhich  results  in  a  slightly  lower  pressure  at  the  boom  inlet; 
settling  time  is  also  shoxt-er,  again  a  reflection  of  the  less  severe 
closure . 

The  linear  closure  resulted  in  a  reduction  of  maximum  transient 
pressure  at  both  the  nozzle  and  boom  inlet;  settling  tine  was 
somewhat  less  than  cur/es  D,  ?,  and  G  bat  greater  than  curve  Z,  The 
reduction  in  pressure  was  due  to  the  constant  rate  of  closure  without 
changes  in  slope. 

Curve  r:  yielded  the  lowest  naximum  transient  pressure  of  the 
group.  Settling  time  was  of  the  same  order  as  the  linear  case.  The 

to 


reduction  ir_  pressure  peaks  was  due  to  she  "cnu::ing”  effect  during 
the  final  stare  of  closure.  This  effect  is  exar.ined  further  in  the 
next  section. 

In  Chapter  VI  curve  D  was  selected  as  the  representative  curve 
shape  for  a  quick  disconnect}  the  results  of  this  section  indicate  a 
sinilar  response  to  curves  D,  2,  ?,  and  1;  therefore,  a  general 
analytical  relationship  applicable  to  all  and  hence,  a  quick 
disconnect  can  be  given  as 


K(t)  =  1  -  T 


where  m  <1 


For  the  quick  disconnect  case,  the  linear  and  snubbed  curves 
are  not  applicable  because  of  the  lower  .maxir.un  pressures  generated. 
Valve  Snubbing  Effects,  Valve  snubbing  was  discussed  in 
Chapter  VI  and  is  defined  for  this  study  as  the  reduction  in  the 
rate  of  valve  area  decrease  during  the  final  ten  percent  of  valve 
area  reduction.  The  effect  of  valve  snubbing  was  investigated  by 
running  the  Laboratory  Test  simulation  with  snub'oed  and  unsnubbed 
curves.  Tine  of  closure  and  steady  state  ccnditions  were  the  sane 
as  the  previous  runs.  The  valve  closure  was  given  by  2q  (l)  and 
K(t)  from,  the  modified  closure  curve  H  as  shovrr.  in  Fig  2oa,  Four 
closures  were  evaluated;  an  unsnubbed  cur"v'e,  and  three  with  snubbing 
of  and  11/t,  respectively.  The  snubbed  curve  is  the 

basic  curve  H,  The  effect  of  the  increases  percentages  of  snubbing 
is  shovm  in  Fig  26b,  Increasing  the  snubbing  percentage  decreases 
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tne  T.axiriur.  transient  pressure, 

.\C-13';  l-.odel.  The  effect  of  closure  ::_rve  ar.d  tir.e  of 

closure  on  the  ?IG-135  syster.  described  ir.  i'ropter  V  ■.-.■es  exar.inei  by 
•.trying  the  tine  of  closure  for  several  curves.  The  '/alve  closure  vas 
siven  by  Eg  (l)  with  ll(t;  from  closure  curves  A,  E,  Z,  Z,  F,  G,  H,  and 


the  linear  curve.  Time  cf  closure  was  -.-ariei  from  C,01  c-nd  0,15  3c-c 

to  encompass  the  minimum  closure  tine  discussed  ir.  C.rapter  VI  and  the 

minimum  used  by  Parks  (Ref  10:53) •  The  surre  boot  precharge  was 

50  psig  and  all  other  conditions  were  as  described  in  Chapter  V, 

The  magnitude  of  the  first  peak  at  the  nozzle  was  chosen  as  the 

basis  for  comparison  of  the  results;  the  peak  nu.ntering  is  illustrated 

in  Fig  27.  The  response  of  the  system  subsecuent  to  the  first  peak 

was  basically  unaffected  by  changes  in  closure  ti.me  curve  shape, 

formally  the  most  important  peak  for  design  purposes  is  the  first  peak 

after  valve  closure,  unless  subsequent  peah.s  are  cf  greater  magnitude. 

The  latter  case  is  possible  in  complex  systems  like  the  KC-135  where 

different  concentrations  of  reflected  pressure  '.raves  coalesce  to 

form  peaks  greater  than  previous  peaks.  The  transient  response  in 

Fig  Zf’i  closure  curve  A,  T  =  O.O5  sec;  is  an  example  cf  such  a  case, 

0 

Shown  in  Fig  28  is  the  pressure  of  the  first  peak  at  the  nozzle 
versus  closure  time  for  each  curve. 

For  all  curves  there  is  an  increasin.g  pressure  vrith  decreasing 
closure  time  relationship,  Cur/es  E,  F,  and  3  are  least  affected  by 
changes  in  closure  time.  These  curves  have  closure  rates  that  are  hi.p 
during  the  latter  stages  of  closure  thus  .mosa  of  the  closure  takes 


place  over  a  short  percentage  of  any  closure  time ,  The  slope  of  the 
curve  does  not  reverse,  i,e,,  go  from  posial'.'e  tc  negative  or  vice  versa 


Fis  27.  I’uiubering  of  Transient  Press’are  Peaks  for  a  Typical 

Response  to  Valve  Closure  for  the  KC-135  '-cciel  using 

curve  A  and  T  =  0,05  sec, 
c 

over  the  closure  interval. 

Curves  3,  C,  and  H  have  similar  responses  to  decreasing  closure 
tine.  All  of  these  curves  appear  to  "be  very  dependent  on  changes  to 
closure  tine,  espscially  below  0,03  seconds,  The  extreme  effects  of 
closure  tine  \'ariation  is  shown  by  the  apprcxirately  20C;-:  increase  i 
pressure  for  a  reduction  in  closure  time  from  0,03  to  0,01  sec. 

The  linear  curve  exhibits  a  relationship  of  pressure  increase 
versus  closure  tine  that  is  similar  to  curves  H,  E,  and  G  except  tha 

the  increase  is  not  quite  so  dramatic.  This  is  a  consequence  of  the 
constant  slope  with  no  snubbing. 
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Fig  2';.  Effect  of  Closure  Tine  Variation  on  1—  Peak 
Transient  Pressure  for  the  KC-135  Model 


-.siHi"  Curve  results  in  a  noclest  cnar.-'e  over  tho  range  of 
closure  ti'cc,  Z:.tz  is  tc  be  expected  cince  the  curve  Is  heavily 
snulfced  and  has  no  reversal  of  Elone,  liote  that  neither  curve  A  nor 


curves  a, 


and  G  have  slope  reversals;  the  effect  of  chances  in 


closure  tine  is  siailar  for  all.  The  lov.er  values  of  curve  A  are 
due  tc  cr.u:hing  el'feots, 

Cverall,  for  t;.e  KG-i35f  effect  of  reduciruj  closure  tLv.e  is 
to  shift  the  dominant  pressure  peak  fron  the  third  peak  to  the  first 
peak  for  all  cuirves  except  Z,  F,  and  G,  First  peak  pressure  for  all 
cases  increases  with  reduction  in  closure  tine. 

The  consistency  of  the  third  peak  value  can  be  attributed  to  the 
system  configuration,  since  for  cor.plex  systems  the  compression  viave 
initiated  by  the  valve  closure  does  not  bring  the  flow  to  rest  as  it 
might  in  simple  systems;  therefore  the  time  required  for  flow  stoppage 
is  characteristic  of  the  system.  As  shewn  by  To  (3) ,  a  ijressure  peak 
will  occur  at  the  point  of  complete  flow  steppage  and  dependent  on  the 
system,  it  is  likely  to  maintain  a  consistent  magnitude  without  regard 
to  closure  time.  For  the  KC-135  model,  complete  flow  stoppage  occurs 
coincident  with  the  timing  of  the  third  peak.  This  is  illustrated  in 
Fia  29  fer  a  typical  response  upstream  of  the  nozzle. 

Parks  (?.sf  10:53“55)  investigated  pressure  versus  closure  time  fo 
the  interval  0,05  --a  0,15  sec  and  found  little  change.  For  the 
closure  used  by  Parks;  curve  A,  =  0,05  sec;  the  results  previously 
discussed  indicate  that  the  third  peak  would  be  predominant  and 
consistently  the  sane  magnitude.  Therefore,  the  results  Parks 


obtained  are  consistent  with  the  results  of  this  study. 


'ip  25.  Typical  Effect  of  rlovr  Stoppa,-^  or.  Trr.r.sient  Responce 

fcr  t'io  KC-135  Ilodel  usin^  curve  A  and  T  ■  0,05  sec 

c 


The  effect  oP  r>ur?'e  fcoot  preche.r7e  inv-sti-r.ted  by  •.'aryin;7 
the  precharre  fi'oi;;  50  '^•0  100  psic*  A  lir.ear  closure  cura'^e  was  used 
with  closure  times  of  0,03  O.O5  sec.  All  other  conditions  were 
as  described  in  Section  V, 

Results,  shov.’n  in  Fig  30,  indicate  ar.  increase  in  pressure  for 
the  first  peak  and  a  reduction  in  pressure  for  the  third  peak  over 
the  range  of  increasing  precharge  pressure ,  Nozzle  pressure  versus 
precharge  pressure  curve  slopes  for  the  first  peak  are  similar  for 
both  closure  tines  with  a  higher  transient  pressure  for  the  0,03  sec 
closure.  The  third  peak  has  the  same  response  for  both  closure  times. 

Parks'  study  (Ref  10:53)  indicated  a  decreasing  maximum 
transient  pressure  with  increasing  precharge  pressure  relationship. 

As  discussed  earlier,  due  to  the  slow  closure  (T^  =  0,05  sec)  and 
extreme  snubbing  of  the  closure  curve  usei  by  Parks,  the  third  peak 
was  predominant  in  his  study.  Therefore,  the  results  obtained  by 
Parks  and  the  third  peak  results  of  this  stuiy  are  in  agreement. 

The  differing  effects  on  the  first  and  third  peaks  result  from 
the  inability  of  the  accumulator  to  respond  guickly  to  short  duration 
impulses.  The  third  peak  pressure  builds  at  a  slower  rate  as  the  flow 
decreases  to  zero.  For  the  first  peak,  increasing  the  precharge 
pressure  acts  to  harden  the  system,  previling  progressively  less 
energy  abserptien  for  the  impulse.  Conversely,  the  slower  building 
pressure  of  the  third  peak  is  attenuated  •■.•inh  the  increased  pre charge 
since  more  energy  is  required  to  compress  the  gas  in  the  accumulators. 
The  effect  of  precharge  for  the  KG-135  -odel  is  consistent  with 
the  results  obtained  for  the  Laboratozp’  Test  riodel.  The  Laboratory 


VIII 


Conclusions 


Icised  on  the  reisearch  and  simulations  conducted  and  reported 
herein,  it  is  concluded  that: 

1 ,  The  HYTfJA.1'  computer  proexam  can  predict  the  transient 
response  of  simple  systems  uith  good  r-'ults, 

2,  oinulation  of  the  Laboratory  Test  resulted  in  a  favorable 
agreement  v;ith  experimental  findings  for  the  general  transient 
response  but  underpredicted  maximum  transient  pressure  by 
approximately  157^. 

3,  For  the  Laboratory  Test  simulation: 

3-)  A  757-  decrease  in  accumulator  volume  vilth  other 
variables  held  constant  had  a  negligible  effect  on  maximum  transient 
pressure, 

b)  Increasing  the  accumulator  precharge  pressure  with  other 
variables  held  constant  caused  an  increase  in  transient  pressure, 

c)  Decreasing  the  diameter  or  increasing  the  length  of  the 
accumulator  entry  lines  with  all  other  variables  held  constant  caused 
an  increase  in  transient  pressure  and  a  decrease  in  the  time  required 
to  settle  to  a  steady  pressure  after  valve  closure, 

4,  The  representative  closure  curve  fcr  quick  disconnects  must 
be  derived  from  analysis  of  the  combined  nozzle/receptacle  dynamics. 
The  curve  developed  in  this  study  gave  good  results  when  used  in  the 
Laboratory  Test  simulations  and  compared  to  the  experimental  results. 

5,  Increased  percentages  of  valve  snubbing  decreases  transient 
pressures;  this  is  consistent  with  the  findings  of  Kinzig  (Ref  8). 


Addition5--l  verification  of  oonputer  simulations  vlth 
exp-srimentally  measured  data  is  needed  tc  =ain  oonfider.ce  in  the 
ccr.r-uter  predictions.  A  scaled  down  system  similar  in  ccnfiruiration 
to  the  Lacoratory  Test  could  be  used  in  the  laboratory  for  parr-neter 
variation  exneriments. 


A  fuel  systems  oriented  deri'/ative  of  hl’TRA!.'  is  available;  this 


program,  FUZLTRAN,  uses  the  basic  HYT3AN  prosram  with  additional 
subroutines  added  for  fuel  system  component  models  (Ref  I5).  FUSLTRAl 


was  not  used  for  this  study  because  confidence  could  not  be  developed 
in  the  results  obtained  from  it,  FUiLTRAh’  should  be  evaluated 
thoroughly  to  determine  what,  if  any,  usage  limitations  exist. 
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